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Abstract 

We calculate the effective Coulomb repulsion between electrons/holes, Um\ and site energy for an 
isolated BEDT-TTF [bis(etliylenedithio)tetrathiafulvalene] molecule in vacuo. U^^ =4.2 ±0.1 eV 
for 44 experimental geometries taken from a broad range of conformations, polymorphs, anions, 
temperatures, and pressures (the quoted 'error' is one standard deviation). Hence we conclude 

(v) 

that Um is essentially the same for all of the compounds studied. This shows that the strong 
(hydrostatic and chemical) pressure dependence observed in the phase diagrams of the BEDT-TTF 

(v) 

salts is not due to Um ■ Therefore, if the Hubbard model is sufficient to describe the phase diagram 
of the BEDT-TTF salts, there must be significant pressure dependence on the intramolecular terms 
in the Hamiltonian and/or the reduction of the Hubbard U due to the interaction of the molecule 
with the polarisable crystal environment. The renormalised value of Um'^ is significantly smaller 
than the bare value of the Coulomb integral: Fq = 5.2 it 0.1 eV across the same set of geometries, 
emphasising the importance of using the renormalised value of Um\ The site energy (for holes), 

(v) 

Cm = 5.0 lb 0.2 eV, varies only a little more than Um across the same set of geometries. However, 
we argue that this variation in the site energy plays a key role in understanding the role of disorder 
in ET salts. We explain the differences between the (3l and (3h phases of (BEDT-TTF)2l3 on the 
basis of calculations of the effects of disorder. 
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I. INTRODUCTION 



It is well known that density- functional theory (DFT), as implemented using the cur- 
rent generation of approximate exchange- correlation functionals, fails to capture the physics 
of strongly correlated electrons, for example, the Mott insulating state.- The electrons are 
strongly-correlated in the layered organic charge transfer salts of the form (ET)2X, where 
ET is bis(ethylenedithio)tetrathiafulvalene (also known as BEDT-TTF and shown in Fig. [1]) 
and X is a monovalent anion. A great deal of interesting physics is driven by these strong 
electronic correlations including unconventional superconductor, Mott insulator, charge or- 
dered insulator, bad metal, charge ordered metal, and pseudogap phases.-i^ To date, first 
principles atomistic theory has not been able to describe the full range of physical phenom- 
ena observed in the ET salts. Therefore most theoretical approaches to describing organic 
charge transfer salts have been based on effective low-energy Hamiltonians, such as Hubbard 
models.-i^ However, in order to make a detailed comparison between these calculations and 
experiment accurate parameterisations of the relevant effective low-energy Hamiltonians are 
required. 

The strong electronic correlation effects evident in the macroscopic behaviour of the ET 
salts do not arise because the intramolecular Coulomb interactions are unusually strong. 
Rather, the electrons move in narrow bands and thus the 'kinetic' energy associated with 
electrons hopping between molecules is much smaller the equivalent energy scale in elemental 
metals; thus the narrow bands give rise to the strong correlations.-i^ This means that density 
functional theory (DFT) is likely to give accurate results for single ET molecules, even if 
it cannot accurately calculate the properties of the extended system. Similar statements 
hold for the alkali doped fuUerides,^ oligo-acene and thiopenes,- and the organic conductor 
tetrathiafulvalene-tetracyanoquinodimethane (TTF-TCNQ)^ where significant insight has 
been gained from parameterising effective low-energy Hamiltonians from DFT calculations 
on a single (or a few) molecules. There have been some pioneering efforts^'^'^'^'^ to apply 
this approach to ET salts, but no systematic study of the wide range of ET salts known 
to synthetic chemistry. Many other families of organic charge transfer salts based on other 
donor and, even, acceptor molecules are also known to show strong correlation effects.— 

In this paper we present DFT calculations for a single ET molecule in vacuum. By 
studying various charge states we calculate parameters for the Hubbard model, specifically 
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FIG. 1: The eclipsed and staggered conformations of the ET molecule. 



the effective Coulomb repulsion between two electrons/holes, 'the Hubbard U\ and the site 
energy. In sectionlTTlwe define the problem and stress the difference between the renormalised 
and bare values of the parameters we calculate below. In section HVl we compare and contrast 
various methods of calculating the Hubbard parameters and provide benchmarks for our 
calculations against the previous literature. In section |V] we report our main results: the 
values of the parameters calculated at 44 distinct geometries observed by x-ray scattering in 
a broad range of conformations, polymorphs, anions, temperatures, and pressures. We find 
no significant variation in the Hubbard U, but a slightly larger variation in the site energy. 
We discuss the importance of this variation of the site energy for the effect of disorder on 
the superconducting state of the ET salts. In particular we are able to explain the origin of 
the difference between the superconducting critical temperatures of the Ph and Pl phase of 
(ET)2l3 on the basis of these results and the theory of non- magnetic impurity scattering in 
unconventional superconductors. The technical details of our DFT calculations are reported 
in section IIIII and our conclusions are drawn in section IVII We will not attempt to solve the 
Hubbard models relevant to the ET salts in this paper; this is an extremely demanding task 
which has been the subject of numerous previous studies.-i^ 

II. THE HUBBARD MODEL AND DENSITY FUNCTIONAL THEORY 

It has frequently been argued^'^^ that the physics of ET salts can be understood in 
terms of (various flavours of) Hubbard models. Here two different classes of ET salt need 
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to be distinguished. In the first class, which includes the k, /3, and (3' phases, the ET 
molecules are strongly dimerised. The intra-dimer dynamics are typically integrated out of 
the effective low-energy Hamiltonian, to leave one with a half-filled Hubbard model, where 
each dimer is treated as a site.^ However, if one is to satisfactorily parameterise this effective 
Hamiltonian, one must first understand the intra-dimer dynamics. In the second class of ET 
salt, which includes the a, and 9 phases, the dimerisation is weak or absent and one is 
forced to consider a quarter-filled Hubbard model in which each site corresponds to a single 
monomer .- 

It is convenient, for our purposes, to write the (extended) Hubbard model in the form 

i ij 

where Tii are the terms that depend on the physics of a single ET molecule and Tiij are the 
terms coupling more than one ET molecule. In this paper we assume that 

T~(-i = Eq + emcl^Ci^ + Umcl-^Ci^cl^Cii, (2) 

where c-P annihilates (creates) a particle with spin cr on site i, Eq may be thought of as 
the energy due to the 'core electrons' and the nuclei [a more rigourous definition is given 
below Eq. ([5])], is the site energy for site i, and Um is the effective Coulomb repulsion 
between electrons on site i, whose interpretation we will discuss presently. (In general 
Eq, em, and Um may be different of different sites, but we suppress site labels in order 
to simplify notation. The subscript m serves to remind us that we are dealing with the 
Hubbard parameters for a monomer rather than dimeric parameters.) Note that we do 
not consider vibronic interactions; however these can straightforwardly be included via the 
Hubbard-Holstein model. Calculations of the vibronic interactions have been reported by 
several groups.—*^ Also note that, we will not discuss the form or the parameterisation of 
Hij below. Most previous studies of Hij have been based on the extended Hiickel or tight 
binding approximation, which neglect all two-body terms in TiijM^ 
It is convenient to separate Um into two terms 

Um = U^^-6U^\ (3) 

where Um^ is the value of Um for an isolated ET molecule in vacuum and 5Um^ is the reduc- 
tion of Um when the molecule is placed in the polarisable crystalline environment. 6Um^ has 
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been successfully calculated for alkali doped fuUerenes,— li^*^ oligo-acene and thiopenes,- and 
TTF-TCNQ.- Perhaps the simplest system for which to calculate 5f/m^ is the alkali doped 
fullerides. This system has (approximate) spherical symmetry and a large intermolecular 
spacing relative to the size of the molecule. 5Um^ has been calculated for the alkali doped 
fullerides both by assuming that the crystal is a dielectric medium and by assuming that 
each molecule in the crystal contributes to 5Um^ only via the classical polarisability tensor 
appropriate to that molecule.— li^ii^ii^S These calculations show that 6Um^ can be a significant 
fraction of Um'' , perhaps as much as a half.— However, the calculation of the polarisa- 
tion correction to the U for the ET salts is a much more difficult problem. First, the systems 
have very low symmetry. Second, many ET salts have polymeric counterions, therefore it 
is not possible to approximate the anion layers by a lattice of dipoles. Third, ET molecules 
are separated by distances typical of vr-stacking in aromatic compounds (~3.5-4 A), i.e., 
distances that are quite small compared to the size of the molecule. This invalidates the 
approximations based on well-separated dipoles or an effective mean-field dielectric medium 
that have been successfully applied to the fuUerenes. Therefore, we limit our attention on 
the calculation of Um^ in this paper. Clearly the calculation of 5Um^ is an important prob- 
lem that must be solved before a full parameterisation of the Hubbard model for ET salts 
can be completed. Interestingly, Cano-Cortes et al.^ have recently reported a calculation of 
5f/^^ for TTF-TCNQ, which manifests some of these problems, where the dipole moment 
is associated with individual atoms, rather than entire molecules. 

It is often stated that the Hubbard U is the Coulomb repulsion between two electrons in 
the highest occupied molecular orbital (HOMO) of an ET molecule.— This is, indeed, often 
a conceptually helpful way to think about the Hubbard U . However, if one were to take this 
literally it would be natural to equate uln' with the zeroth Slater-Condon parameter, i.e., 
the Coulomb integral. 



where pa{f) the density of spin a electrons at the position r in the HOMO of the ET 
molecule. Fq is rather different from the Hubbard U . To see this recall that the Hubbard 
model is an effective low-energy Hamiltonian and therefore the Hubbard parameters are 
renormalised. Assuming that Um^ = Fq is equivalent to simply ignoring the high energy 
degrees of freedom not contained in the Hubbard model. This is clearly inadequate: Fq is 
the unrenormalised value of Um\ In principle one should derive the effective low-energy 
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Hamiltonian by explicitly integrating out the high energy degrees of freedom. However, 
this is not practicable in this context for all but the simplest systems.— Therefore, one 
should construct the effective low-energy Hamiltonian so as to capture the relevant degrees 
of freedom and parameterise this effective Hamiltonian appropriately. We now review how 
to do this for the Hubbard model. 

It is trivial to solve the single site Hubbard model ([2]); one finds that 

El = Eo + e^ (5a) 
E2 = Eo + 2em + U^ (5b) 

where En is the energy of the model with n electrons. Note that En may be calculated in 
DFT as it is simply the ground state total energy of ET(2-")+ and DFT does well for total 
energies. It follow immediately from ([5]) that 

Um = Eo + E2 ~ 2Ei. (6) 

This can be visualised as the energy cost of the charge disproportionation reaction 2(ET^) — >■ 
ET + ET^^ for infinitely separated ET molecules. Another interpretation is that Um is the 
charge gap, i.e., the difference between the chemical potentials for particles and holes, for 
ET+. It is also clear from Eq. ([5]) that the site energy is given by 

^m = El — Eq. (7) 

This is simply the second ionisation energy of the ET molecule, ET^ — > ET^^. 

One can also write Eq. ([2]) in terms of the hole operators defined by h]^ = Cia which gives 

'Hi = E2 + ^nihl^hia- + Umh^^^hi^h\Jlil, (8) 

where, 

= -(e™ + Um) =Ei~E2 (9) 

is the site energy for holes. Physically corresponds to the ionisation energy of ET. Thus 
one sees that while the Hubbard U is the same for electrons and holes, the site energy is not. 
Note that we limit our study to salts of the form ET2X where X is a monovalent anion. 
In these materials the HOMO of the ET molecule is, on average, three quarters filled with 
electrons or, equivalently, one quarter filled with holes. Therefore, the hole parameters are 
of more physical relevance than the electron parameters, as we will discuss further below. 

6 



One might also consider the Taylor series for the dependence of the energy on the 'clas- 
sical' (i.e., not quantised) charge in terms of electrons, 



E{q) ^ E{0) + — 



, 1 d^E 



q=2- 



or holes, 



E{q) ^ E{2) + — 



9=0+ 



2 dq" 



1 d^E 
^^2 5^ 



(2 - qf 



q=2' 



(10a) 



(10b) 



q=0+ 



where q = 2 indicates the limit as g — 2 from below and q = indicates the limit as 
q from above. Again, E{q) can be calculated from DFT by calculating the total energy 



of the relevant charge states. It then follows from 

d^E 



that the classical U is given by 



dq^ 



Similarly, the classical site energies are 



dE 
dq 



and irnc = -TT- 

q=2- dq 



(11) 
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q=0+ 



as one expects from Janak's theorem.— We, therefore, are left with the classical results in 
terms of electrons, 

E{q) ~ E{0) + e^e(2 - q) + ^Umc{2 - qf. (13a) 

and holes, 

E{q) ~ E{2) + UcQ + \u^,q\ (13b) 

Note, however, that = — (cmc + '^Umc) in contrast to Eq. ([9]). This is a manifestation of 
the failure of classical approximation for the site energy. 

The density functional formalism can be applied to fractional numbers of electrons.-*^ 
Therefore one can calculate the energy for a large number of charge states and find the best 
fit values of -E(O), e^c, ^mc, and U^c in Eqs. f fT3|) . which is operationally more satisfactory 
than determining the two free parameters (^^ and Um determine e^) from two data points 
[E{1) — -E(O) and E{2) — -E(O)]. Of course this method suffers from the fact that, unlike 
Eqs. (E]), (CD, and (El), Eqs. ^ are not exact. 

Eqs. ( |T3l) give remarkably good fits to the fractional charge data (c.f. Fig. [2]). Further 
we find that Uml is an extremely accurate approximation to Um (typically the two values 
differ by < 1%). In contrast e^c and are not good approximations to and This 
is a consequence of our use of an approximate generalised gradient approximation (GGA)^- 




FIG. 2: (Color online) Dependence of the energy of an ET molecule on the charge of the 
molecule. Two different geometries are studied. The 'frozen' geometry is relaxed from that found 
experimentally^ in /3-(ET)2l3 in the neutral charge state and then held fixed during the self- 
consistent field (SCF) calculations at different charge states. For the 'relaxed' data the nuclear 
geometry is relaxed separately for each charge state. The energies are exactly equal by definition 
in the charge neutral state; for other charge states the nuclear relaxation lowers energies. This 
reduces the curvature of E{q), i.e., lowers Uml and hence Um\ It can be seen that for both the 
'frozen' and the 'relaxed' geometries the energy of the fractional charge states is extremely well 
described by the classical quadratic functions, Eqs. (jl3p . In contrast it is known that for the 
exact functional the energy of molecules with fractional charges is a linear interpolation between 
the integer charge states.— '^^'^^ We find that to an excellent approximation Umc=Um^ whereas emc 
and ^mc are not good approximations to em and ^m- In contrast for the exact functional em = ^mc, 

(v) 

Cm = £,mc, and Umc = 0. This incorrect result is a manifestation of the delocalisation error of DFT 
and closely related to the band gap problem.— i^^i^I These calculations use TZP basis sets and TM2 
pseudopotentials. 

exchange-correlation functional. It is known that for the exact functional the energies of 
states with fractional charges are linear interpolations between integer charge states.-'^^i'^ 
In contrast for local functionals, such as the local density approximation (LDA) or that of 
Perdew, Burke and Erzenhof (PBE),^^ E{q) is a convex function. ^'^^ This is closely related 
to the self-interaction error and leads to the delocalisation error in local functionals and 
the band gap problem .->2^>2I Indeed, for the exact functional e.^ = emc, Cm = Cmc, and 
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III. COMPUTATIONAL METHODS 



We used DFT to calculate the total energies and HOMO charge densities of ET molecules 
in various conformations and charge states. In order to compare the effects of different basis 
sets and pseudopotentials (including no pseudopotential) we made use of two different im- 
plementations of DFT, SIESTA-- and NRLMOL.— In each case a separate self- consistent 
field (SCF) calculation was performed at every charge state studied. 

SIESTA calculations were performed using the PBE exchange-correlation functional.— A 
triple-C plus single polarisation (TZP) basis set of Sankey numerical atomic orbitals^^ was 
used for all atoms. These orbitals are confined to some radius from their centres, which 
introduces a small increase to the energy of each orbital. The value of was determined 
by specifying the maximum allowed increase in energy due to this cutoff, which we limited 
to 2 mRy. The convergence of the integration mesh is determined by specifying an effective 
plane-wave energy cutoff for the mesh, which we set to 250 Ry. The initial spin moments 
on each atom were arranged antiferromagnetically, i.e., with opposite signs on neighbouring 
atoms, wherever possible. All SIESTA calculations reported below used pseudopotentials 
constructed according to the improved TrouUier- Mart ins (TM2) method.— 

NRLMOL performs massively parallel electronic structure calculation using Gaussian 
orbital methods. We again used the PBE exchange correlation functional. In the calculations 
presented in this paper we use Porezag-Pederson (PP) basis sets^'^ which have been carefully 
optimized for the PBE-GGA energy functional using variational energy criteria. As discussed 
in Ref. 



30l . for each atom, the basis sets are optimized with respect to the total number 
of Gaussian decay parameters, and with respect to variation of these parameters and the 
contraction coefficients. These are roughly of triple to quadruple zeta quality.— As compared 
to other Gaussian basis sets, a key improvement in the PP optimization scheme is that the 
resulting basis sets satisfy what is now referred to as the Z^^^^ theorem. This theorem^S 
discusses proper scaling of the Gaussian exponents near the nuclei as a function of atomic 
charge. It has been shown that the resulting PP basis sets exhibit no superposition error and 
alleviate the need for counterpoise corrections in weakly bound systems. We report both all 
electron calculations and calculations using pseudopotentials for C and S atoms performed in 
NRLMOL. These pseudopotentials were constructed using the method of Bachelet, Hamann 
and Schliiter (BHS).^ 
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Nuclear positions for C and S atoms were obtained from x-ray 
crystallography.— i^i^i^i^SdSilid^ilSJiJ^J&dlil^d^i^i^ In section |IV] we will discuss vari- 
ous schemes for relaxing these atomic positions and for calculating Um\ therefore we delay 
defining these methods to that section. The total energy was then computed for charge 
states ranging from charge neutral (doubly-occupied HOMO), to the 2+ charge state (the 
same orbital unoccupied) in increments of 0.1 electrons. The resulting 21 data points per 
molecule were then fitted to a quadratic function by the method of least squares, from 
which values of Uml, emc and were extracted according to Eqs. (fT3l) . Um-, Cm and 
were calculated from integer charge states according to Eqs. (EI), (171), and ([9]). 

Calculations of the Coulomb integral, Fq, make use of the orbitals from the corresponding 
calculations. The integrals in Eq. (jl]) were calculated in the charge neutral state on a 
55 X 30 X 20 mesh using the trapezoidal rule. We have also investigated finer meshes and 
the calculations were found to be well converged with respect to the number of integration 
points. Equation (jll) has a pole at ri = Y2- Corresponding terms in the numerical integral 
were approximated by the analytical solution of equation (jl]) for a sphere of uniform charge 
density and the same volume as the mesh volume elements. These points contribute ~ 2% to 
the value of Fq, and the error associated with the uniform sphere approximation is estimated 
to be < 1% of the contribution from individual mesh points, i.e. < 0.02% of the total value 



IV. GEOMETRY OPTIMISATION, PSEUDOPOTENTIALS, BASIS SETS, AND 
CALCULATION METHODS 

We report results of calculations of and Fq performed using TZP basis sets and TM2 
pseudopotentials in Table [H For these calculations we take the nuclear positions of the C 
and S atoms from x-ray scattering experiments^ on /5-ET2I3 and the nuclear positions of 
the H atoms were relaxed using the conjugate- gradient method.— Henceforth, we will refer 
to this method of determining the nuclear geometry as 'experimental'. Both eclipsed and 
staggered conformations (c.f. Fig. [1]) are observed experimentally^ and we present results 
for both conformations. However, we find no significant differences between the eclipsed and 
staggered conformations in either f/m'* or Fq. We also find no significant difference in the 
calculated values of Um^ and Uml ■ Similar results are found at the other geometries studied 
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Quantity 


Eclipsed 


Staggered 


rji-) 


4.12 


4.07 




4.17 


4.08 




5.40 


5.35 



TABLE 1: Comparison of different methods of calculating the effective monomer on-site Coulomb 

(v) (v) 

repulsion for ET. No significant differences are found between Umc and Um from the DFT cal- 
culations. This result has been confirmed for all of the geometries studied below. Further, no 
significant differences are found between the two conformations. Fq is significantly larger than 
Um\ in qualitative agreement with previous results from wavefunction based methods (c.f. Table 
In]). The nuclear geometry is that seen 'experimentally' in /3-(ET)2l3 for both the eclipsed and 
staggered conformations. The calculations use TZP basis sets and TM2 pseudopotentials. All 
values are in eV. 

below. 

Several groups^i^'^'i^ have reported calculated values of Fq or Um^ previously. We sum- 
marise their results in Table [Tll It can be seen from Tables [T] and |TT] that there is broad 
agreement between the values calculated by a range of different levels of theory, basis sets, 
and starting from different experimental geometries. Across all of the calculations Fo is 
~40% larger than Um/' ■ This over-estimation has been attributed to the fact that Fq ne- 
glects orbital relaxation,-i^i^ i.e., the tendency of the orbital to change as its population 
changes. Conceptually, this is very important. It stresses that the sites in the Hubbard 
model do not correspond to any particular 'orbital' of the ET molecule, as these orbitals 
change in response to charge fluctuations, etc. in the crystal. Fq is, in fact, simply the first 
order approximation to Um^ Therefore the true value of Um^ contains many non-trivial 
quantum many-body effects which are absent in Fo beyond just orbital relaxation.— 

Now we turn to investigate the effect of different calculation schemes on the Hubbard 
parameters. To that end, results are presented for values of Um^ for ET molecules in geome- 
tries taken from /5-(ET)2l3, calculated with two major modifications to the calculation. The 
first modification is to relax the nuclear coordinates. We have investigated three methods of 
determining the nuclear coordinates.— The first nuclear co-ordinate set is the 'experimental' 
geometry defined above, i.e., the experimental positions are used for the C and S atoms, 
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Method 


Basis set 


Atomic geometry from 






SBK-31G* 


K-(ET)2Cu[N(CN)2]Cl^ 


5.90 


HFl 


SBK-31G* 


K-(ET)2Cu[N(CN)2]Br69 


5.90 


HF^. 


SBK-31G* 


K-(ET)2Cu[N(CN)2]I^ 


5.83 


RHF {q = 4}^ 


6-31G** 


K-(ET)2Cu[N(CN)2]Bria 


5.44 


RHF {q = 2}^ 


6-31G** 


K-(ET)2Cu[N(CN)2]BrIS 


6.40 


RHF {q = 0)8 


6-31G** 


K-(ET)2Cu[N(CN)2]Br22 


5.58 


Mean Fq 






5.84 


Std. dev. 






0.33 


Method 


Basis set 


Atomic geometry from 




RHF^ 


6-31G** 


K-(ET)2Cu[N(CN)2]Br™ 


4.48 


VB-HF^ 


AMI 


K-(ET)2Cu(NCS)2^^ 


3.92 


VB-HF^ 


AMI 


K-(ET)2Cu(NCS)2^ 


3.83 


VB-RHFi° 


AMI 


/3-(ET)2l3^ 


3.90 


VB-RHF^o 


DZ 4-31G 


/3-(ET)2l3^ 


4.61 


Mean U^^ 






4.15 


Std. dev. 






0.37 



TABLE 11: Previously reported calculations of the Coulomb energy (Fo and U^') of ET. These 
calculations were performed at various levels of theory, using various basis sets and with geometries 
taken from x-ray crystallography experiments on various different materials. Fq is significantly 

(v) 

larger than Um , consistent with our results. We make use of the following abbreviations in the 
above table: HE (Hartree-Fock) ; RHF (restricted Hartree-Fock) ; andVB (valence bond). Fortunelli 
and Painelli® calculated Fq for ET dimers in different charge states, q. All energies are in eV. 

but the H atoms, which are not seen in x-ray scattering experiments, are relaxed. The 
second set was obtained by relaxing all the nuclear positions in the charge neutral state and 
holding that geometry fixed for all later SCF calculations at different charge states. This 
set is labelled 'frozen'. The third set of nuclear co-ordinates was obtained by relaxing all 
nuclear positions in each charge state. This set is labelled 'relaxed'. The second modifi- 
cation to the computational method is to perform the calculations with different basis sets 



12 



and pseudopotentials. 

We present results for the experimental, frozen, and relaxed geometries in Table UTTl using 
a variety of basis sets and pseudopotentials. The TZP basis sets consistently gives ~ 0.2 
eV smaller than the PP basis sets; similar small differences are seen in and However, 
the trends between the different geometries are reproducible between the two calculations. 
These differences may be related to the fact that the PP orbitals are optimised for each 
charge state whereas the TZP numerical orbitals are the same for every charge state. This 
additional degree of freedom reduces the curvature of E{q) which reduces Uml and hence 
Um^ . Whether pseudopotentials are used or all electron calculations are performed does not 
have a significant effect on the results. The experimentally measured^ gas phase ionisation 
energy of ET is 6.2 eV, in reasonable agreement with our calculated value of C,m- Note, 
however, that ^mc is significantly smaller than the experimental value. 

The trend between the different geometries is more interesting. Each time a geometry 
relaxation is performed Um^ is reduced. This results form the increased number of degrees 
of freedom when more geometry relaxation is allowed. This can be seen from the plot 
of the energy for different charge states in Fig. [2l For both of these curves both the 
geometry relaxation and SCF calculation were performed using the TZP basis sets and 
TM2 pseudopotentials. The two q = results are therefore, in fact, the same data point. 
As the charge is varied away from q = the geometry remains fixed in 'frozen' calculation 
but relaxes in the 'relaxed' calculation, thus the energy of 'relaxed' data is lower than that 
of the 'frozen' calculation g > 0. This effect gets larger as q increases, thus Um^ is smaller 
in the 'relaxed' data than for the 'frozen' calculation. The physical content of this result 
is simply the fact that intramolecular vibronic couplings act to lower Um^ However, one 
should expect this effect to be rather stronger in vacuo than in the crystalline environment 
where the motion of the molecule is significantly constrained. 

Thus we must face the question: should one use the 'experimental', 'frozen', or 'relaxed' 
geometry to calculate f/m'*? We believe that the 'experimental' geometry gives the most 
useful information. Firstly, there are small differences in the reported geometries for different 
ET salts, and one would like to understand the effect of these on a single molecule before 
considering the effects of changes in the crystal structure on the emergent physics of the 
crystal. Secondly, the experiments effectively 'integrate over' all of the relevant charge 
states and therefore provide an 'average' conformation. Thirdly, the experiments naturally 
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Coordinates 


Conformation 


Basis set 


Pseudo. 


^ m 






^ mc 






RMSE 


experimental 


staggered 


PP 


All 


3.90 


-9.15 


5.25 


3.87 


-11.07 


3.33 


0.010 


experimental 


staggered 


PP 


BHS 


3.93 


-9.19 


5.26 


3.90 


-11.13 


3.34 


0.012 


experimental 


staggered 


TZP 


TM2 


4.08 


-9.28 


5.21 


4.07 


-11.31 


3.17 


0.012 


frozen 


staggered 


PP 


All 


3.78 


-9.63 


5.85 


3.76 


-11.50 


3.99 


0.009 


frozen 


staggered 


PP 


BHS 


3.81 


-9.68 


5.87 


3.78 


-11.57 


4.00 


0.010 


frozen 


staggered 


TZP 


TM2 


4.11 


-9.28 


5.17 


4.09 


-11.33 


3.14 


0.010 


relaxed 


staggered 


TZP 


TM2 


4.07 


-9.07 


5.00 


4.00 


-11.08 


3.08 


0.019 


Mean 


staggered 






3.95 


-9.33 


5.37 


3.92 


-11.28 


3.44 




Std. dev. 


staggered 






0.13 


0.24 


0.34 


0.13 


0.20 


0.39 




experimental 


eclipsed 


PP 


All 


3.99 


-9.40 


5.40 


3.97 


-11.38 


3.45 


0.012 


experimental 


eclipsed 


PP 


BHS 


4.03 


-9.46 


5.42 


4.00 


-11.45 


3.45 


0.013 


experimental 


eclipsed 


TZP 


TM2 


4.17 


-9.26 


5.09 


4.12 


-11.32 


3.08 


0.013 


frozen 


eclipsed 


PP 


All 


3.78 


-9.63 


5.85 


3.76 


-11.50 


3.99 


0.009 


frozen 


eclipsed 


PP 


BHS 


3.81 


-9.68 


5.87 


3.78 


-11.57 


4.00 


0.010 


frozen 


eclipsed 


TZP 


TM2 


4.10 


-9.33 


5.23 


4.07 


-11.36 


3.22 


0.011 


relaxed 


eclipsed 


TZP 


TM2 


3.98 


-9.04 


5.05 


3.95 


-11.01 


3.12 


0.010 


Mean 


eclipsed 






3.98 


-9.40 


5.42 


3.95 


-11.37 


3.47 




Std. dev. 


eclipsed 






0.14 


0.22 


0.33 


0.14 


0.18 


0.38 





TABLE III: Calculated bare and renormalised parameters for the Hubbard model for ET monomers 
under various geometry relaxation schemes and with different pseudopotentials, basis sets and codes 
(see section |HI|) . The 'experimental' geometry is that reported for an ET molecule in /3-(ET)2l3r^ 
measured at 298 K, with the H atom (not observed in x-ray crystallography) positions relaxed. The 
'frozen' coordinate system was relaxed in the charge neutral state and held fixed for other charge 
states. The 'relaxed' geometry was optimised at every charge state. All geometry relaxations 
were carried out in calculations using TZP basis functions and TM2 pseudopotentials; we have 
also carried out the relaxations in using the other methods in the table and find no significant 
differences. The abbreviation pseudo. (for pseudopotential) is used in this table and others below. 
The RMSE is taken from the fit to the classical Eqs (|13|) . All values are in eV. 
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include the effects on the molecular conformation due to the crystalline environment, which 
are absent from in vacuo calculations. Therefore we now move on to consider the value of 
Um^ found for the 'experimental' geometry found in x-ray scattering experiments from a 
wide range of different ET salts. 



V. IMPURITY SCATTERING, POLYMORPHISM AND CHEMICAL PRES- 
SURE 



We begin by studying polymorphism in (ET)2l3, for which the experimental literature 
contains more reported polymorphs than any other salt of ET. In Table HVl we report the cal- 
culated values of the Hubbard parameters and the Coulomb integral for ET molecules in 'ex- 
perimental' geometrie82^i^>^i^i2£ taken from the a, (3, n, and 6 polymorphs of (ET)2l3. All 
of the calculations in this section were performed using TZP basis sets and TM2 pseudopo- 
tentials. There is little variation in Um^ across any of the polymorphs {U^^ = 4.22±0.09 eV). 
For the a, k, and 9 polymorphs is extremely uniform (^^ = 4.840 ± 0.008 eV) but C,m is 
somewhat larger in the (3 polymorph. More interestingly is 0.12 eV larger for the stag- 
gered conformation of /?-(ET)2l3 than it is for eclipsed conformation (the two conformations 
are sketched in Fig. [1]). To understand the significance of this result we must first briefly 
review a few experiments on /?-ET2l3 and a little of the theory of impurity scattering in 
unconventional superconductors. 

The application of hydrostatic pressure, P, to (3-ET2I3 has a dramatic effect on the 
superconducting critical temperature, Tc. At ambient pressure ~ 1.5 K but when the 
applied pressure reaches P ~ 1 kbar a discontinuous increase in (~7 K) is observed. The 
low Tc state (P < 1 kbar) is labelled the /5l phase and the high Tc state (P > 1 kbar) is 
labelled the /5h phase. When the pressure on the /5h phase is decreased the material does 
not return to the /9l phase but rather Tc is seen to further increase. Below T ~ 130 K 
the resistivity of the /3h phase is found to undergo a discontinuous decrease while no such 
anomaly is found in the /5l phase. Incommensurate lattice fluctuations have been observed 
in the /3l phase but they are absent in the jSn phase below T ~ 130 K.— The incommensurate 
lattice fluctuations are stabilised by variations in the conformational ordering of the terminal 
ethylene groups of the ET molecules and thus can only exist in the presence of disorder.— 
For a more detailed review of this phenomenology see Ref. [12 
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Crystal type 


conformation 


Fo 


rjiv) 






^ mc 






RMSE 




staggered 


5.15 


4.27 


-9.10 


4.83 


4.19 


-11.22 


2.84 


0.027 




staggered 


5.14 


4.29 


-9.13 


4.84 


4.20 


-11.25 


2.85 


0.030 




staggered 

OO ------ 


5.35 


4.08 


-9.28 


5.21 


4.07 


-11.31 


3.17 


0.013 




eclipsed 


5.40 


4.17 


-9.26 


5.09 


4.12 


-11.32 


3.08 


0.013 




eclipsed 


5.22 


4.31 


-9.16 


4.85 


4.22 


-11.29 


2.85 


0.030 




eclipsed 


5.40 


4.21 


-9.05 


4.84 


4.15 


-11.14 


2.84 


0.022 


Mean 




5.28 


4.22 


-9.16 


4.94 


4.16 


-11.26 


2.94 




Std. dev. 




0.12 


0.09 


0.09 


0.16 


0.05 


0.07 


0.15 





TABLE IV: Calculated bare and renormalised parameters for the Hubbard model for an ET 
molecule at the 'experimental' geometry observed at ambient temperature and pressure in var- 
ious polymorphs of (ET)2l3. The changes in conformation due to the crystal packing structure do 

(v) 

not have a large effect on the value of Um ■ The calculated is larger for the (3 polymorph than 
the others explored, indeed varies by less than 1% among the other polymorphs. Note that 
there is a significant difference between the values of in the eclipsed and staggered conforma- 
tions in the /3 phase. This is consistent with the effects of conformational disorder on /3-ET2l3J^ 
The calculations use TZP basis sets and TM2 pseudopotentials. All values are in eV. 

It is well known that non-magnetic disorder can lead to the suppression of in uncon- 
ventional ('non-s-wave') superconductors, as the ET salts are believed to be.-*^ The details 
of this suppression are described by the Abrikosov-Gorkov equation,— 

where T^o is the superconducting transition temperature of a pure sample, 1/r is the quasi- 
particle scattering rate and ip is the digamma function. If we consider only the contributions 
to r from terminal ethylene group disorder then 

-=Ns'kD{Ef)\^U\\ (15) 
r 

where, is the number of molecules in the staggered conformation, D{Ef) is the density 
of states at the Fermi energy, and A^m = 0.12 eV is the difference between the values of 
calculated in the staggered and eclipsed phases.— It has previously been shown^ that 
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over both the (3h and /3l phases of /5-ET2I3 the critical temperature scales with the residual 
resistivity, which is proportional to the contribution to 1/r from impurity scattering, as 
predicted by the Abrikosov-Gorkov equation. 

Our finding that the staggered and eclipsed conformations have a small but important 
difference in C,m shows that the disorder in the conformational degrees of freedom of the 
terminal ethylene groups of the ET molecules, required to stabilise the lattice fluctuations, 
is responsible for the suppression of Tc in the jSi phase. The terminal ethylene disorder 
increases the quasiparticle scattering rate and that this leads to the observed different T^s 



63 



shows 



of the Pu and /3l phases. A simple calculation based on the derivation in Ref. 
that the calculated value of accounts for the observed differences between the critical 
temperatures and residual resistivities of the Ph and Pl phases with a few percent of the 
molecules staggered conformation. 

It is interesting that such signiflcant differences are caused by the conformation changes 
in ET. In Fig. [3] we plot the Kohn-Sham orbitals corresponding to the HOMOs of eclipsed 
and staggered ET molecules at the 'experimental' geometries found in /3-ET2I3. Note that 
there is very little electron density on the terminal ethylene groups, which are the only parts 
of the molecule in different positions in the 'experimental' eclipsed and staggered phases. 
Furthermore, there are no signiflcant changes between the HOMO electron densities on the 
other atoms between the two conformations. 

It is well known that changing the anion in ET salts (which changes the unit cell volume) 
has a remarkably similar effect to applying a hydrostatic pressure.- Therefore, changing the 
anion is often referred to as 'chemical pressure'.-*^ Both hydrostatic and chemical pressure 
have dramatic effects on the phase diagram of ET salts. This is typically understood in terms 
of a variation of Hubbard model parameters, e.g., the ratio U/t, where t is the relevant 
hopping integral.-!^ Therefore, we now consider the changes in the Hubbard parameters 
caused by changing the anion. Figs. HI [H and [6] show the calculated values of, respectively, 
Um\ €m, and C,m foT 'experimental' geometries taken from x-ray scattering experiments on 
crystals with a wide range of anions and several different polymorphs. The crystallographic 
data sets^^i^'^i^i^i^i^i^i^i^'^'^'^'^i^'^'^ used to determine the nuclear positions include 
data taken at a range of temperatures and under different pressures. Details of the calculated 
Hubbard parameters, and the crystallographic measurements on which the calculations are 
based are given in the supplementary information. 
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(a) 






FIG. 3: (Color online) Orthographic projections of the isosurface of the HOMO of the ET molecule 
in 'experimental' geometries taken from (a) eclipsed /3-(ET)2l3, (b) staggered /3-(ET)2l3, and (c) k- 
(ET)2Cu[N(CN)2]Cl. Note the great similarity of the HOMOs corresponding to the 'experimental' 
geometries from eclipsed and staggered /3-(ET)2l3 and the small electronic density on the terminal 
ethylene group, which is involved the change between the eclipsed and staggered conformations. 
Further, the HOMOs of the ET molecule in the 'experimental' geometries taken from /3-ET2I3 
and K-(ET)2Cu[N(CN)2]Cl are remarkably similar despite the fact that this geometry is taken 
from a different crystal polymorph with a different anion. This is consistent with our finding that 
the changes in the conformation of the ET molecule, in different polymorphs and in crystals with 

(v) 

different anions, do not significantly affect Um ■ Colour indicates the sign of the Kohn-Sham orbital. 
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All isosurfaces are ±0.07 A~^/^ and calculated in the charge-neutral state. These calculations use 
TZP basis sets and TM2 DseudoDotentials. Animations showine different isosurfaces are available 



Anion 
A AuBr^ 
□ Aulj 

■ PC0I3 

▼ Cu[N(CNyBr 

V Cu[N(CNyCI 

G Cu[N(CNyi 

• Cu(NCS)2 

+ I3 

• lAuBr 
X IBrj 

A (SF),(CH),(CFyC0)3 

• (ET),Hg,3,Br3 



FIG. 4: (Color online) The effective Coulomb repulsion between electrons/holes, Um' , on an 
ET monomer in the 'experimental' geometries corresponding to different anions, conformations 
(eclipsed offset to the left, staggered to the right), temperatures, pressures, and crystal poly- 
morphs. Note the limited range (3.9 - 4.4 eV) of the ordinate. We see that Um^ does not change 
significantly across the different (ET)2^ crystals and has a mean value of 4.2 it 0.1 eV. The cal- 
culations use TZP basis sets and TM2 pseudopotentials. Full details of the parameterisation are 
given in the supplementary information. 



It can be seen from Fig. 0] that the none of the changes to the crystal (anion, polymor- 
phism, conformation, temperature, or hydrostatic pressure) have a significant effect on Um^ , 
indeed we find that for all of the structures we have studied Um^ = 4.2 ± 0.1 eV (here, and 
below, the reported 'error' is simply one standard deviation in the data and does not re- 
flect systematic errors, particularly those arising from our approximate, density functional, 
calculation of the total molecular energies). 

An important question is: how large does variation of U/t have to be in order to explain 
the experimental results? We can estimate this from strongly correlated theories of the 
Hubbard model. Both the resonating valence bond theory,— and cluster extensions to dy- 
namical mean- field theory^ suggest that U/t is required to vary by perhaps a factor of two 
in order to explain the experimental phase diagram, although variational quantum Monte 
Carlo results suggest that the variation need not be quite this large.— 

Our results show that f/m^ is a transferable property and does not play any significant 
role in the 'chemical pressure' effects observed in the behaviour of the ET salts. This means 
that either the chemical pressure comes from 6Um^ and/or the elements of Tiij, or else the 





▲ 
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K 
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Crystal polymorph 
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Anion 
A AuBr^ 
□ Aulj 

■ PC0I3 

▼ Cu[N(CNyBr 

V Cu[N(CNyCI 

G Cu[N(CNyi 

• Cu(NCS)2 

+ I3 

• lAuBr 
X IBrj 

A (SF),(CH),(CFyC0)3 

• (ET),Hg,3,Br3 



FIG. 5: (Color online) The site energy for electrons, Cm, on an ET monomer in the 'experimental' 
geometries corresponding to different anions, conformations (eclipsed offset to the left, staggered to 
the right), temperatures, pressures, and crystal polymorphs. The mean value is = — 9.2ib0.1 eV. 
The calculations use TZP basis sets and TM2 pseudopotentials. Full details of the parameterisation 
are given in the supplementary information. 



simple Hubbard model description is not sufficient to describe the behaviour of the ET salts. 

Figs. [5] and [6] show that the variation of em = —9.2 ±0.1 eV and = 5.0 ± 0.2 eV is 
not much larger than that in Um'' . However, it is known that the ET salts are very sensitive 
to disorder—*^ therefore these subtle effects may lead to much more profound effects in the 
extended system than the variation in Um\ Note also that while the Hubbard U is the same 
for electrons and holes the site energies are rather different, as expected from the discussion 
in section [Tll when U is large. 

The similarity between the electronic states of ET molecules in the conformations found 
in different salts is further emphasised by fact that the Kohn-Sham orbital corresponding to 
the HOMO is essentially the same in all of the molecules we have studied. As an example of 
this we plot the electronic density of the HOMO for an ET molecule in the 'experimental' 
geometry taken from a crystal of k-ET2Cu[N(CN)2]C1 in Fig. [3l This can be seen to be 
remarkably similar to the equivalent plots for /3-ET2I3 (also shown in Fig. [3]) despite the 
fact the polymorph and the anion are different. 




Crystal polymorph 
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Anion 
▲ AuBr^ 
□ Aulj 

" PC0I3 

▼ Cu[N(CNyBr 

V Cu[N(CNyCI 

O Cu[N(CNyi 

• Cu(NCS)2 

+ I3 

♦ lAuBr 
X IBrj 

A (SF),(CH),(CFyC0)3 

Crystal polymorph * (^"'")4'^92.89^'^8 

FIG. 6: (Color online) The site energy for holes, ^m, on an ET monomer in the 'experimental' ge- 
ometries corresponding to different anions, conformations (eclipsed offset to the left, staggered to 
the right), temperatures, pressures, and crystal polymorphs. shows only slightly more variation 

(v) 

than Um ■ However, small changes in are known to have significant effects on the supercon- 
ducting state observed in ET salts.— 1^ (Because these ET molecules are quarter filled with holes 
in the salt, rather than em, is the relevant site energy to consider when discussing the role of 
disorder.) The mean value is = 5.0 it 0.2 eV. The calculations use TZP basis sets and TM2 
pseudopotentials. Full details of the parameterisation are given in the suppmentary information. 



VI. CONCLUSIONS 



In summary, we have calculated the Hubbard U and site energy for an isolated ET 
molecule. We found that Um^ = 4.2 ±0.1 eV for 'experimental' geometries taken from a 
broad range of conformations, polymorphs, anions, temperatures, and pressures. That is, 
Um^ is essentially the same for all of the compounds studied. 

The dependence of the macroscopic behaviour of the ET salts on (hydrostatic and chem- 
ical) pressure is usually understood in terms of the variation of the Hubbard parameters 
with pressure. Strongly correlated theories of the Hubbard model^i^i^ suggest that U/t is 
require to vary by a factor of ~2 in order to explain the experimentally observed chemical 
pressure effect. Therefore the pressure dependence must be contained in either the cor- 
rection to Um from the crystalline environment, 6Um\ or the intermolecular terms in the 
Hamiltonian, Ti^. The renormalised value of Um^ is significantly smaller than the bare value 
of the Coulomb integral, Fq = 5.2 ±0.1 eV across the same set of geometries, emphasising 
the importance of using the renormalised value of Um\ 
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The site energy (for holes) = 5.0 ± 0.2 eV varies only a little more than Um across 
the same set of geometries. However, we have argued that this variation plays a key role 
in understanding the role of disorder in ET salts in general and in explaining the difference 
between the Pl and Ph phases of /5-ET2I3 in particular. 
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